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SUMMARY 


The method developed in NASA TN D-319 for studying the atmosphere 
entry of vehicles with varying aerodynamic forces has been applied to 
obtain a closed-form solution for the motion, heating, range, and varia- 
tion of the vehicle parameter m/CpA for nonlifting entries during which 
the rate of increase of deceleration is limited. The solution is appli- 
cable to vehicles of arbitrary weight, size, and shape, and to arbitrary 
atmospheres. Results have been obtained for entries into the earth’s 
atmosphere at escape velocity during which the maximum deceleration and 
the rate at which deceleration increases were limited. A comparison of 
these results with those of NASA TN D-319, in which only the maximum 
deceleration was limited, indicates that for a given corridor depth, 
liEiiin§ the rate of increase of deceleration and the maximum deceleration 
requires an increase in the magnitude of the change in m/CpA and results 
in increases in maximum heating rate, total heat absorbed at the stagnation 
point , and range . 


INTRODUCTION 


Man’s performance of useful duties during atmosphere entry requires 
that the maximum deceleration and the rate at which the deceleration 
increases be held within prescribed limits (see refs. 1, 2, and 3). The 
problem of limiting the maximum deceleration has been studied by several 
investigators. In the study reported in NASA TN D-319 (ref. 4 ), a dif- 
ferential equation applicable to vehicles of arbitrary weight, size, and 
shape, and to arbitrary atmospheres was developed for entries during 
which the aerodynamic forces vary. A closed-form solution for nonlifting 
entries was obtained for the motion, heating, range, and variation in 
m/CpA required to maintain specified maximum deceleration limits. The 
purpose of the present investigation is to obtain a solution of the 
equation developed in reference 4 for the case of specified deceleration- 
rate limits as well as specified maximum decelerations. The solution is 
used to determine the effect of limiting the deceleration and the rate of 
increase of deceleration (by varying m/CpA) on the corridor depth, 
heating characteristics, and range for vehicles making shallow nonlifting 
entries into the earth’s atmosphere at escape velocity. 
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NOTATION 


reference area for drag, sq ft 

2D 


drag coefficient, 


pV 2 A 


.-2 


-2 


drag force, lb 

gravitational acceleration, ft sec 
gravitational conversion constant, 32*2 ft sec 
deceleration in g units 

specified rate of increase of deceleration, •—•, g/sec 

CL U 

lift force, lb 

mass of vehicle, slugs 

convective heating rate per unit area at the stagnation point, 
Btu ft " 2 sec -1 

dimensionless function proportional to convective heating rate 
(eq. (26)) 

total convective heat absorbed per unit area at the stagnation 
point, Btu ft " 2 

dimensionless function proportional to total heat absorbed, 
(eq. (28)) 

distance from planet center, ft 

radius of curvature of vehicle surface, ft 

range measured from entry point, ft (see eq* (29)) 

time, sec 

tangential velocity component normal to a radius vector, 
ft sec " 1 


dimensionless velocity ratio, 


resultant velocity, 


u 


cos 7 


dimensionless velocity ratio, 


u 

■s/S? 


V 

v/gr 


ro o v n 
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w 

weight of vehicle at earth T s surface, mg G , lb 


y 

altitude, ft 


z 

dimensionless function of u (eq. (k)) 


3 

atmosphere density decay parameter, ft" 1 , 1/23 500 ft" 
earth 

1 for 

7 

flight -path angle relative to the local horizontal, negative 
for descent 

A 

vehicle parameter, — — 

g d a 


Ay p 

corridor depth, statute miles 


p 

atmosphere density, slugs ft -3 


Po 

mean value for exponential approximation to atmosphere 
altitude relation, slugs ft" 3 , 0.0027 slugs ft" 3 for 

density- 

earth 


Subscripts 


i 

initial 


max 

maximum 


mod 

modulated entry 


0 

surface of a planet 


unmod 

unmodulated entry 


i, begin 

point in the trajectory vhere modulation begins 


2 

point in the trajectory where modulation for constant decelera- 
tion rate ends 

3 , end 

point in the trajectory where modulation ends 


I, II, III 

phase of a trajectory 



Superscript 

differentiation with respect to u 
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ANALYSIS 


In the present analysis of shallow nonlifting entries, m/CpA will 
be varied so that both the maximum deceleration and the rate at which the 
deceleration increases do not exceed specified limiting values. The 
trajectory for this type of modulated entry generally consists of the 
four distinct phases illustrated below. The deceleration-time histories 
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for an unm odulated entry and a modulated entry are shown for the same 
initial conditions. During phase I the vehicle enters the atmosphere with 
a given velocity (ui) and entry angle (7i) and maintains a constant 
m/CpA. The deceleration increases monotonically until the specified value 
of the deceleration rate is attained at velocity u x , at which velocity the 
magnitude of the deceleration is Gi. At this velocity phase II begins 
with the modulation of m/CpA in such a manner that the deceleration 
rate is held constant at a specified value K; that is. 


G 


= Gi 



■fci) 


for t x < t < t 2 


(l) 


ir\ o cm 



5 


The velocity u 2 , at which phase III begins, is determined as that 
velocity at which the specified value of the maximum deceleration G 2 
is first attained. During phase III, m/Cj)A is modulated in such a 
manner that the deceleration is held constant until the velocity is reduced 
to the value u 3 * This velocity (u 3 ) is determined by specifying that 
termination of the modulation will allow the deceleration to decrease 
monot oni c ally as the vehicle completes its entry. Phase IV is this final 
phase of the entry. It should be noted that phases I, III, and IV of this 
report correspond to phases I, II, and III, respectively, of reference k. 

The general differential equation which defines the trajectory for 
the four phases of the modulated entry (and the unmodulated entry) is 
obtained from reference k as 


uZ Tt 




+ uZ 


All 

A 



+ 




cos 4 / 


0 


(2) 


where 


is the dimensionless independent variable. 



2A 


(3) 


( 4 ) 


is the dimensionless dependent variable first presented in reference 5 
and 


A s 


m 

C D A 


A(u) 


The accompanying expression for the flight -path angle is 


(5) 


Vp r sin 7 = Z ' - % + Z ( 6 ) 

Equation (2) is applicable to vehicles of arbitrary weight, size, and 
shape, and to arbitrary atmospheres. The solution of equation (2) for the 
present analysis consists primarily of determining the solution for 
phase II of the modulated portion of the trajectory, since the solution 
for phase III is indicated in reference 4 and the solutions for the unmod- 
ulated portions of the trajectory (phases I and IV) are indicated in 
reference 5 * 


There are four basic steps to the solution of equation (2) for 
phase II. First, by eliminating Z and its derivatives from equation (2) 
there results a differential equation which defines the required 
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variation in m/CpA. Second, the solution of this differential equation 
is obtained. Third, the duration of phase II Is determined (i.e., the 
velocity u 2 is determined) . Finally, the so Lution for phase II is 
matched with the solutions for phases I and Ii::. 

The first step is accomplished by means of the equation from reference 
4 for the resultant deceleration during shallow nonlifting entries 1 
which is 


G = Sr uZ (7) 

An expression for the Z function during phase II can be obtained by 
eliminating G from equation ( 7 ) in the following manner. The first 
derivative of equation (l) can be written 


— = K = — — 

dt du dt 


(8) 


Substitution of dG/du, as obtained from differentiation of equation ( 7 ), 
and du/dt, as obtained from reference A for l/d = 0, 

= -s/Pg uZ (9) 

in equation (8) yields the differential equation for the Z function 
during phase II 


A 


uZ(uZ' + Z) = - 


If it is noted from equation (7) that 


K 


*'/pg •Vf 


(10) 


t r 


ai£ 


= 2 uZ(uZ> -- Z) 

du du 


( 11 ) 


and that at u » Ui, (uZ) 2 = equation ( 10 ) can be integrated to 

yield the desired expression for Z during phf.se II 


Z = \/a + bu 

u 


:i2) 


where 


x In this report a shallow entry implies that during modulation the 
flight-path angle is sufficiently small that cos y ~ 1 is a valid 
approximation . 


in o oj 
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a = 1 - bui 


b 


2K 

Vp g Vpt" b x 2 


Bi 


Gi 

n/Pt 


> (13) 


Substituting equation (12) and its derivatives in equation (2) gives the 
differential equation which describes the variation in m/cpA required 
to maintain a constant deceleration rate 


Af 

A 



b 

_2(a + bu) 



+ 


k_ 

-2 

u 


— --3 b . bf If Ui! . . = 0 (Ik) 

2u( a + bu) k(a. + bu) 2 B x 2 (a + bu) 


and the related expression for the flight-path angle is obtained from 
equation (6) as 


sin 7 - 


Bi -s/a + bu 
s/Pr 


u 


L2(a + bu) u 


_ 2 A^_ 

~ “ A 


(15) 


The solution of equation (lU) for the variation of m/Cj)A, or A , 
with velocity during phase II can be obtained in closed form as 

$Xi ■ exp ksW x(s)(2a - frliT - 


^aC 
3b 2 L 


X(g) - 1 


2Cui 

3b 


|X(g) - 1 


ka 


,b 2 Bi 2 


+ 2 ) Zn | + 


gui 

3bBi £ 


(1 - I) 


2au 2 


v (1 - t z ) + (1 - £ 3 ) 


3b 2 Bi 2 ' 3 ' 9bB x 2 


2(3a - l) Zn V5 + 1 


3 b 2 B! 2 Jg - 


(16) 



u 


vnere 




x(g) = Ji - bu a (i - i) 


> (17) 


C .2E sin n + 

■Dl 


Zn 


•fa + 1 2(3a - l) 


/aBi 2 yfei - 1 


3b 2 B x 2 


The details of the derivation of equation (l6) are presented in the 
appendix to this report. The expression for the flight-path angle obtained 
from equation (15) and the derivative of equation (l6) is 


s ±n 7 IX = ~ 


Bi 

■TPr 


2n 


■Ja + X(g) _ 4a 


sfa Bi 2 s[a - X( | ) 3b 2 Bi S 


X(S) + 


2u, 


3bB x " 


SX(g) - C 


(13) 


The velocity u 2 at which phase II is terminated is determined from 
the deceleration-velocity relation as that value of velocity at which the 
specified maximum deceleration G 2 is first attained. From equations 
( 7 ), (12), (13), and (17), the deceleration-velocity relation is found to 

be 


A 


G = Gi Jl - bui(l - gf (19) 

and from equation (19), the value of g 2 (and hence u 2 ) at which phase II 
is terminated is 



For modulated entries during which only the deceleration rate is limited, 
the velocity u 2 is determined as indicated ir the appendix. 

It now remains to match the solution for jhase II with the solutions 
for phases I and II. This is done by matching the attitude and flight - 
path_angle at the velocities for the beginning and end of phase II (ux 
and u 2 , respectively). The quantities at the beginning of modulation 
Zi, or A 1? and are known from phase I (obtained as indicated in 
ref.J?); hence, with the aid of equation (15) and the fact that at 
u = a + bu = 1, the solutions for phases I and II are matched by 
satisfying the conditions 


ro o vji 



A IIi = A 1 
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•JOr — 

— u x sin y ± 



b 

2 


( 21 ) 


where the subscript Up denotes initial values for phase II. At the 
ena of phase II, A 2 , and y 2 are known (eqs. (lb) and (l 8 )), and the Z 
function and the equation for the flight-path angle are known for phase III 
(see ref. ), hence, the solutions for phases II and III are matched by 
satisfying the conditions 


where 


A IIIi = A 2 



— - u 2 sin 7 
B 2 


_ 2 _ 

u 2 


( 22 ) 


B 2 


^/pr 


(23) 


For modulated entries during which only the deceleration rate is limited 
the solutions for phases II and III are matched as indicated in the 
appendix . 


The solution of equation ( 2 ) is known for the remainder of the 
trajectory; that Is, the closed-form solution for the variation of m/CpA 
required to maintain constant maximum deceleration during phase III is 
given in reference 4 and the solution for phase IV is given in reference 5 . 
However, It should be noted that in the present analysis, for phase III 



where (A 2 /a 1 ) ii is determined from equation ( 16 ) with £ = £ } and 
(A/AsJjjj is determined from the closed-form solution of reference 4 . 
As noted earlier, the closed-form solution of reference 4 for the 
variation of A required to maintain a constant deceleration is for 
phase II of that reference. 


RESULTS AND DISCUSSION 


The closed-form solution of this report and that of reference 4 have 
been employed to study the effect on nonlifting entry trajectories of 
maintaining a specified deceleration rate and a specified maximum 
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deceleration by varying W/CpA during entry into the earth's atmosphere 
at escape velocity. 2 In what follows, the effect of modulation on the 
time history of W/CpA, flight-path angle, deceleration, and velocity is 
discussed first. The effect of modulation on the corridor depth, heating, 
and range is then discussed. 


Trajectory Parameters 


Typical time histories of W/CpA, flight-path angle, deceleration, 
and velocity for nonlifting vehicles are presented in figure 1 for unmod- A 

ulated and modulated shallow entries into the earth's atmosphere at 5 

escape velocity (Vi = 1.4) . The curves shown are for an entry angle 0 

7 i = -6 . The modulated entries are for a lOg maximum deceleration limit 2 

only (for which (dG/dt) max = O.95 when G = 10 ), and for the combined 


limits of a 0.25g/sec deceleration rate and a lOg maximum deceleration. 

It is interesting to note that the maximum deceleration rate for an unmod- 
ulated lOg entry (y ± = -4.32°) is 0.22g/sec. Tice data in figure 1 for 
the unmodulated entry and the entry modulated to limit the maximum decel- 
eration only are reproduced from reference 4. For the entry considered 
above, it can be seen in figure l(a) that addin;; the restriction of a 
0.25g/ sec deceleration-rate limit to a lOg maxljium deceleration limit 
requires that the W/CpA change capability of vhe vehicle must be 32 
rather than 21. By beginning the modulation at an earlier time to satisfy 
the added restriction of the deceleration-rate limit, the flight -path 
curvature is slightly reduced relative to that for the entry which 
satisfies only the maximum deceleration limit (see fig. 1 (b)) and this 
results in the reduced deceleration rate shown In figure l(c) during 
phase II of the modulation period. 


Corridor Depth 


In this report the overshoot boundary Is defined by the trajectory 
for which a vehicle would pass through just enough atmosphere to reduce 
the velocity to local circular as the vehicle is. about to exit from the 
atmosphere (see ref. 5 ). The undershoot boundai-y is defined by the 
trajectory for which a specified deceleration rate and/or maximum 
deceleration is not exceeded. For vehicles for which the initial w/OpA 
Is identical for both the overshoot and undershoot boundaries, the 
corridor depth, for successful entry into an exponential atmosphere, 
depends only upon the allowable entry angle (see ref. 4). Changes in 
corridor depth considered herein result from the use of drag modulation 
to maintain specified deceleration time historic s during entry along the 

2 Since the remainder of the report is concerned only with entries 
into the earth's atmosphere W/CpA rather than m/CpA is used. Also, 
for earth, the value of 30 is used for -v/^r, and 27 is used for l/ -/pg. 
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undershoot boundary. For example, increased corridor depths result from 
lowering the undershoot boundary by the use of modulation which permits 
entries at steeper angles than are possible without modulation. 

Corridor depths and the corresponding entry angles attainable by 
modulating W/CpA during shallow nonlifting entries into the earth's 
atmosphere at escape velocity are shown in figure 2 as a function of the 
ratio of W/CpA at the end of modulation to w/CpA at the beginning of 
modulation. Results are shown for modulated entries for which the maximum 
deceleration is limited to 5.0, 7.5, and 10. Og. For these entries the 
maximum deceleration rate is not limited and increases with increases in 
entry angle. For the above specified deceleration limits, results are 
shown which indicate the effect of imposing the additional restrictions 
of deceleration -rate limits of 0.10, 0.25, and 0.50g/sec. With the 
restriction of a deceleration-rate limit of O.lOg/sec the maximum decel- 
eration for entry angles less than -5.O3 0 is less than lOg, and the maxi- 
mum deceleration for entry angles less than -4.32° is less than 7«5g. 
Therefore, the data in figure 2 (and in all subsequent figures) represent- 
ing these conditions are faired with a short dashed curve. It can be seen 
in the figure that the additional restriction of a deceleration-rate limit 
generally reduces the corridor depth. For example, a vehicle with a 
change capability in W/CpA of 21 can enter the earth's atmosphere at 
escape velocity through a corridor up to 30 miles in depth ( 7 ^ - -6°) and 
not exceed a maximum deceleration of 10. Og. However, the additional 
restriction of a deceleration-rate limit of 0.25g/sec reduces the corridor 
depth to approximately 27 miles, which corresponds to an entry angle of 
approximately -5.8°. On the other hand, in order to attain a 30-mile 
corridor and comply with this additional restriction, the change capability 
in W/CpA must be increased from 21 to 32. It is interesting to note 
from reference 6 that to use lift to attain a lO.Og-limited, 30-mile 
corridor, without modulation, the entering vehicle must have an l/d 
capability of at least ±0.3. For this unmodulated entry along the under- 
shoot boundary, however, the deceleration rate reaches a maximum value of 
0.42g/sec. To limit an unmodulated lifting entry along the undershoot 
boundary to 10. Og and 0.25g/sec, the L/D must be restricted to 0.05. 

These restrictions (-O.3 < l/d < 0.05, K max = 0.25, and G mflY = 10), with 
no modulation, provide only a 17-mile entry corridor. This entry corridor, 
however, could be increased by modulating L/D. 


Heating 


For the entries considered herein, peak heating occurs during the 
high-altitude portion of the trajectory where the density and the attendant 
Reynolds numbers are low; consequently, the following discussion of heating 
during entry will cover only those results obtained when a laminar 
boundary layer is considered. 
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Convective heating rate at the stagnation point ,- The convective 
heating rate per unit area at the stagnation point during entry into the 
earth* s atmosphere, as given in reference 5; is 


% 

where for entries considered here with limited values of deceleration 
rate and maximum deceleration 


= 590 nr g 

-Je cJcifiR cos3 r 


q. = 


u 5/2 zi' 2 for phases I and IV 
B 1 1/2 u 2 ( a + bu) 1/4 for phase II 
B 2 1 ' / ‘"u 2 for phase III 


(26) 


Equation (25) has been used to calculate the stagnation -point heat -transfer 
rate for modulated and unmodulated entries which have the same maximum 
deceleration. Results for the maximum heating rate as a function of 
corridor depth are shown in figure 3 for deceleration limits of 7.5 and 
10. Og, and for deceleration-rate limits of 0.10, 0.25, and 0.50 g/sec. 

Also shown are the values of maximum deceleration rate ( dG /dt ) max for 
the entries which are modulated to satisfy only a maximum deceleration 
limit. In making these heating-rate calculaticns it was assumed that the 
radius of curvature R of the vehicle surface at the stagnation point 
was identical and constant for all entries. It was also assumed that the 
Initial w/CpA was the same for all entries. It should be noted that 
the data presented in figure 3 can be used to calculate the ratio of 
stagnation -point heating rates for the same modulated and unmodulated 
entries (same for which R mo d ^ Runmod* provided R remains constant 

during entry. This is done by multiplying the ratio of heating rates 
shown in figure 3 by ( R unm od/^ mo d) 1 2 • If R varies during entry, the 
stagnation-point heating rate will depend upon the variation of R with 
time or velocity (see eqs. (25) and (26)). 


The heating -rate penalty incurred by the use of modulation to increase 
corridor depth can be seen in figure 3* For example, the maximum heating 
rate for the lOg -limited modulated entry (lOg only) which provides a 
30-mile corridor is about twice that for the unmodulated entry which pro- 
vides a 7 _n iile corridor. Adding the further restriction of limiting the 
deceleration rate to 0.25g/sec results in a maximum heating rate of almost 
two-and-a-half times that for the unmodulated entry. 

Total convective heat absorbed at the stagnation point .- The total 
convective heat absorbed per unit area at the stagnation point during a 
modulated entry into the earth's atmosphere can be written (see ref. 5 ) 


ro o vji > 
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where 


= 


_ 15,900 




w 


‘ CpAR \du. 


^du 


(27) 


u 


3/2 


Z 1/2 cos 2 7 


for phases I and IV 


aa = y 

du ' 


u 2 


•D 17 2 / \l/4 

Bi (a + bu) 


for phase II 


) ( 28 ) 


—2 

U 

n 1/2 


for phase III 


Equation ( 27 ) has been used to calculate the total heat absorbed per unit 
area at the vehicle stagnation point for modulated and unmodulated entries 
which have the same maximum deceleration limit* The results as a function 
of corridor depth are shown in figure 4 for deceleration limits of 7*5 
and 10. Og, and for deceleration-rate limits of 0.10, 0*25, and 0 * 50 g/sec. 
In all calculations the same assumptions were made with regard to vehicle 
shape and V/CpA that were made in calculating the heating-rate data of 
figure 3 (i.e., R is Identical and constant during entry, and the initial 
W/CpA Is Identical for all entries) . For vehicles which satisfy these 
assumptions, the penalties in total heat absorbed at the stagnation point 
can be seen in figure 4. For example, in the case cited previously, In 
which modulation is used to increase the depth of a lOg-limited corridor 
from 7 to 3^ miles, the total heat absorbed at the stagnation point during 
the modulated entry which limits only the maximum deceleration is 1.8 
times that for the unmodulated entry. Adding the further restriction of 
limiting the deceleration rate to 0.25g/sec results in the absorption of 
2*3 times as much total heat as that for the unmodulated entry. As in 
the case of the heat -transfer rate, the data of figure 4 may be used, 
similarly, to calculate the total heat absorbed at the stagnation point 
for different vehicles where R mo d = constant 7 ^ Runmod; otherwise, R(u) 
must be available for use in equation ( 27 )* 


As noted above, the comparison of the heating characteristics of 
modulated and unmodulated entries was based on the assumption that the 
initial w/GpA was the same for all entries. If, instead, it Is assumed 
that the final w/cpA is the same for all entries, the heating per unit 
area at the stagnation point will be less for the modulated entries than 
for the unmodulated entries. Specifically, the data presented in figures 3 
and 4 will be reduced approximately by the appropriate ratio of 
(^begin/Asnd) 1 2 (i.e., the ordinates of figures 3 and 4 will be less than 
or equal to unity) . The total heat absorbed by the vehicle, however, 
obtained by Integrating the local heat absorbed over the vehicle surface 
may be greater or less for the modulated entries than for the unmodulated 
entries, depending upon the vehicle geometry during entry. 
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Range 


In reference 5 the range between any tvo points of the entry 
trajectory is defined as the cireuinf erential iistance traveled between a 
point n where the dimensionless velocity is u n and a point n + 1 
where it is reduced to Thus, from reference 5 





1 r Un cos 7 


du 


( 29 ) 


In the present analysis equation (29) is used to determine the range 
during phases I and IV. For cos 7 - 1 , Z(u) as given by equation (12), 
and G(u) as given by equation ( 19 ), equation ( 29 ) can be evaluated to 
yield the range during phase II as 


( 30 ) 

Similarly, for Z(u) as given by equation ( 7 ) with constant G, the range 
during phase III Is 


(£5) 

— a 

_ u 2 

[ 1 . 

Ait’ 

Vo/ in 

2G^ 


W . 


Equations ( 29 ), (30), and. ( 31 ) have been used to calculate the range for 
the same modulated and unmodulated entries for -which the corridor depth 
and heating characteristics have been discussed previously. The results 
as a function of corridor depth are shown in figure 5 for G max limits 
of 7*5 and 10. Og, and for the added restrictions of K max limits of 0.10, 
0*25, and 0.50g/sec. It can be seen that the range is generally reduced 
as a result of using modulation to increase corridor depth. However, as 
shown in figure 5(b) for a lOg-limited corridor, imposing the low decelera- 
tion rate K max = 0.10 actually increases the range relative to that for 
the unmodulated lOg-limited entry. 


As 

r o. 


II 


3b 2 Gx 2 


(2a - bu^Gg - (2a - bux)Gx 


A 


CONCLUDING REMARKS 


The differential equation developed in N/lSA TN D -319 for entries 
during which the aerodynamic coefficients and vehicle shape vary has been 
used to obtain a closed-form solution for the motion, heating, range, and 
variation in W/C-qA for nonlifting entries during which the rate of 
increase of deceleration is limited. The solution is applicable to vehi- 
cles of arbitrary weight, size, and shape, and to arbitrary atmospheres. 
The solution was used to calculate trajectory parameters for shallow 
entries into the earth's atmosphere at escape velocity. As an example of 


tn O CM 
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the results obtained, the depth of a lOg-limited corridor can be increased 
from 7 miles for an unmodulated entry to 30 miles for a modulated entry 
for a vehicle with a w/cpA change of 21 if no consideration is given to 
limiting the rate of increase of deceleration. However, if in addition 
to limiting the maximum deceleration to lOg, the maximum rate of increase 
of deceleration is limited to 0.25g/sec, a W/CpA change of 32 is 
required to obtain the 30 -mile corridor* For a given corridor depth, 
limiting the rate of increase of deceleration and the maximum deceleration 
results in increases in maximum heating rate, total heat absorbed at the 
stagnation point, and range, compared with limiting only the maximum 
deceleration* The heating results were obtained for the case in which 
the radii of curvature at the stagnation point were constant and equal 
during all entries, and the initial values of w/opA were identical. 


Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., May 1, 196 1 



APPENDIX 


CLOSED-FORM SOLUTION FOR m/CjjA., u 2 , AND r FOR MODULATED SHALLOW 
NONLIFTING ENTRIES FOR WHICH LG/ It IS CONSTANT 


The variation in m/UpA with velocity required to maintain a constant 
rate at which the deceleration increases Is defined by equation (l 4 ). 

Upon substitution of the dependent variable 


^1 = 1 dA ^ 1 aU/CpA) 

A A du m/cpA du 


equation (l 4 ) reduces to 


_b 1 p a i ■ u 2 


2 (a + bu) uj Bi 2 (a 4 bu) 4 (a f bu) 2 2u(a 4 bu) ^ 


The corresponding expression for the flight-pith angle is (see eqs . (Al) 
and (15)) 

nr: \/a T bu f , b 2] 


7 “ . tg va + bu , b 2 

pr sin 7 ~ Bx — — — p 4 — — — = 

u L 2(1 + bu) u_ 


Integration of equation (A2) yields ( \1 a + buAi is an integrating factor) 


A 1 ^ ^ - 1 u 2 n " n/ a + ^b ; I + 4a 

v a Bi 2 %/a 4 bu Va 4 \/a + bu 3b 2 Bi 


4 a — ki — 2 

5 2 U 2 U 

Bi 2 3 bBx 2 


+ b u + | b +c 

2a (a 4 bu) u 2a 


4 bu 


where C is the integration constant which c in be evaluated in terms of 
the velocity and flight -path angle at the beginning of modulation; that 
is, at u = u 1; a 4 bu = 1 , and from equation (A3), 


P i = Sln/l " 2 + ^ 
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Thus from equation (A4) evaluated at u = u x 


C = 


v3r 


1 


Zn 


•Va + 1 


^ sin + 

Bi i # B x 2 ^ - 1 


Proceeding by integration of equation (Al) 


2(3a - 1) 
3b 2 B x 2 


(A5) 


^ ~~ C^A ~ ex B(*^*P ^u + constant) (A6) 

Equation (A4) is readily integrable and the integration constant can be 
evaluated in terms of the velocity and m/CpA at the beginning of modu- 
lation (u x and A x , respectively) . Thus, the closed-form solution for the 
■variation of m/CpA which maintains a constant rate at which the decel- 
eration increases during a shallow nonlifting entry is 

A _ m/ CpA f o r~ w * v 

g ° ’ exp x(s)(2a - b ^ )in f ffuj - !nXU) 


4aC 

3b 2 


X(s) - 1 


2Cu x 

3b 


£X(§) - l 


+ 


4a 

-3b 2 B x 2 



Zn | 


where 


2ui 

3bB x ' 


'1 - I) - 


2 au x c 

3b 2 B, 


(1 - I 2 ) + 


2u x 3 

9bB x £ 


(1 - t 3 ) 


2 (3a ~ l) , s/a + 1 1 
3 \/ab 2 B x 2 n ^-l J 


t 


u 


x(|) = 7i - bu x (i - 4 ) 

and a, b, and are given by equation ( 13 ). 


(AS) 


The velocity u 2 at vhich modulation should be terminated vhen 
consideration is given to limiting only the rate at vhich the deceleration 
increases is determined in the following manner. For a given set of 
entry conditions and a specified limiting deceleration rate, the constants 
in equation (A 7 ) are fixed, and it is a simple matter to determine from 
this equation that m/Cj)A increases to a maximum value and then decreases. 
Continuation of the modulation after m/CjyA has reached a maximum vould 
sustain a constant deceleration rate above that vhich vould result if the 
modulation vere terminated. Consequently, the value of the velocity u 2 
is determined as that value for vhich A/a^ is a maximum; that is, the 
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value of velocity obtained vhich satisfies t'ue equation which results 
from equating the right-hand side of equation (A4) to zero. When it is 
desired to terminate this modulation at an earlier velocity at vhich a 
specified maximum deceleration is attained, ns in the present analysis, 
the velocity u 2 is determined as indicated in the ANALYSIS section of 
this report. 


The expression for the flight-path angle during the modulation period 
for constant deceleration rate is obtained from substitution of equation 
(A4) into equation (A3) as 


sin 7 


1 Zn Va + X(|) 
•Jpr LVa Bx 2 <JS - x(|) 


4a 

3 b a B : 2 


X(6) + 


3bBx J 

(A9) 


A 


For entries during vhich it is required that only a specified 
deceleration rate shall not be exceeded, the end of the modulation phase 
is matched vith the final unmodulated phase as follows. The Z function 
and A are knovn at the end of the modulation period for constant decel- 
eration rate (eqs. (12) and (l6)), and A x /A - p = 0 at velocity u 2 ; 
hence, vith the aid of equation (A3) and the fact that the flight -path 
angle for an unmodulated entry is given by (oee ref. 5) 


s/Wr sin 7 - Z 1 - \- 
v. 


(A10) 


the solution for the modulated deceleration-rate phase (phase II) and the 
unmodulated final phase (phase III) are matched by satisfying the condi- 
tions 


'illi = B i 


■Ja + bu 2 
u 2 


Z » 


Illi = " Z III 1 


u 2 2(a + bu 2 )J 


' 

\ (All) 


Equations (All) together with the velocity v 2 constitute the initial 
conditions from which the trajectory characteristics during the remainder 
of the entry are calculated as indicated in inference 5. Where it is 
desired to terminate this modulation at an earlier velocity at which a 
specified maximum deceleration is attained, as in the present analysis, 
the various phases of such a trajectory are matched as indicated in the 
ANALYSIS section of this report. 
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Figure 1.- Continued 
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Figure 1.- Concluded 






Figure 2.- Effect of varying w/CpA on_the deceleration- and deceleration-rate-limited entry 

corridor; V± = 1.4, L/D = 0, earth. 
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Figure 4.- Total lajninar convective heat absorbed per unit area at the stagnation point associated 
with the use of modulation along the undershoot boundary of a deceleration- and deceleration- 
rate-limited entry corridor; Vi = 1.4, L/D = 0, earth. 




Figure 4.- Concluded 
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Figure 5*- Range associated with the use of modulation along the undershoot boundary of 
deceleration- and deceleration-rate-limited entry corridor; V± ~ 1*^ L/D = 0, earth. 




Figure 5.- Concluded 




